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Abstract 

Improved bounds on the copula of a bivariate random vector are computed when partial 
information is available, such as the values of the copula on a given subset of [0, 1]^, or the 
value of a functional of the copula, monotone with respect to the concordance order. These 
results are then used to compute model-free bounds on the prices of two-asset options which 
make use of extra information about the dependence structure, such as the price of another 
two-asset option. 
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1 Introduction 

A (two-dimensional) copula is a function C : [0, 1]^ — s> [0, 1] with the following properties: 

i. Boundary conditions: C(0, u) — C{u, 0) = and C(l, u) = C{u, 1) — u for all u G [0, 1]. 

ii. C is 2-increasing, i.e. for every < mi < U2 < 1 and < wi < W2 < 1, one has 

C{U2,V2) + C(mi, «i) - C{UUV2) - C{U2, Vi) > 0. (1) 

The classical Frechet-HoefFding bounds on the distribution function of a two-dimensional ran- 
dom vector, can be expressed in terms of the copula C of this vector: 

W{u, v) :— max(0, u + v — 1) < C{u, v) < min(M, v) :— M{u, v). (2) 

In the presence of additional information on the dependence between the components of the vector, 
these bounds can be narrowed. Nelsen et al. compute the improved bounds when a measure 
of association such as Kendall's r or Spearman's p is given, and the Bertino's family of copulas 
[I] yields best possible bounds when the values of the copula on the main diagonal are known. 
More generally, given a nonempty set of bivariate copulas 5, Nelsen et al. [T^ introduce pointwise 
best-possible bounds of S: 

= sup{C(u,w)|C e 5} and B{u,v) = ini{C{u,v)\C € S}. 

These bounds are in general not copulas but quasi-copulas, and a fortiori they do not necessarily 
belong to the set S. 

In the theoretical part of this paper (section we first compute the improved Frechet bounds 
when the values of the copula on an arbitrary subset of [0, 1]^ are given, and provide a sufficient 
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condition for each bound to be a copula, and therefore, be the best possible bound. This generalizes 
the findings of 10 on the improved Frechet bounds for copulas with given diagonal sections. 
Next, we compute the best-possible bounds when the value of a real- valued functional of the 
copula, monotone with respect to the concordance order and continuous with respect to pointwise 
convergence of copulas is given, extending the results of |2] . 

Since the work of Rapuch and Roncalli [111 it is known that the prices of most two-asset options, 
when the marginal laws of the two assets are fixed, become monotone functionals of the copula 
with respect to the concordance order. The classical Frechet-Hoeffding bounds therefore lead to 
model-free price estimates for such options [11] |4l |3] . 

In sectional we obtain a new representation for the price of a two-asset option, allowing to use 
a quasi-copula. This representation enables us to compute (in section [5]) the improved model- free 
estimates of the option's value when the prices of all single-asset options on each of the two assets 
are known and some extra information about the dependence structure. This extra information 
may be, for example, the price of a different two-asset option (for example, zero-strike spread 
options are often quoted in the market), or the correlation of two assets. This is similar in spirit to 
a recent work by Kaas et al. [5] who compute worst-case bounds on the Value at Risk of a portfolio 
of two assets when the marginals and a measure of association are known. 

2 Preliminaries 

In this section, we recall several useful definitions and results and fix the notation for the rest of 
the paper. In the definition of quasi-copula [2], the 2-increasing property ([T]) is replaced by weaker 
assumptions: 

Definition 1. A (two-dimensional) quasi-copula is a function Q : [0,1]'^ — > [0,1] with the following 
properties: 

i. Boundary conditions: (5(0, u) — Q{u, 0) = and Q{1, u) — Q{u, 1) — u for all u £ [0, 1]. 

ii. Q is increasing in each argument. 

in. Lipschitz property: \Q{u2,V2) — Q{ui,vi)\ < \u2 — ui\ + \v2~vi\ ioi all {ui,vi,U2,V2) e [0,1]^. 

We denote the set of all copulas on [0, 1]^ by C and the set of all quasi-copulas by Q. The 
concordance order is the order on Q defined by Qi -< Q2 if and only if Qi{u) < Q2{u) Vu G [0, 1]^. 
It is clear that all quasi-copulas satisfy the Frechet-Hoeffding bounds ([2]). Similarly, we say that 
Q" — )■ Q pointwise if Q"{u) Q{u) Vu G [0, 1]^. The Lipschitz property implies that in this case 
the convergence is uniform in u. 

For a copula or a quasi-copula C and a rectangle R — [ui,U2] x [wi,W2] C [0,1]^, we define 

VciR) ■■= Ciu2,V2) + CXui,Vi) - C{UI,V2) - C{U2,V1). 

A subset S C [0, 1]^ is called increasing if for all (ai,6i) G S and (02,62) G S either oi < 02 
and 61 < 62 or ai > 02 and 61 > 62. It is called decreasing if for all (ai,6i) G S and (a2,&2) € S 
either oi < 02 and 61 > 62 or ai > 02 and 61 < 62. It is easy to see that for a decreasing set S, 
the set S := {(a, b) : (a, 1 — &) G S} is increasing. In the same spirit, if C is a copula, the function 
C(m, v) :— u — C{u, 1 — u) is also a copula and if Q is a quasi-copula, Q{u, v) := u — Q{u, 1 — w) is 
also a quasi-copula. 

The following well-known result (see e.g. Theorem 3.2.3 in [5]), gives the best-possible bounds 
of a set of copulas taking a given value at a given point. 

Proposition 1. Let C be a copula and suppose C{a,b) = with (a, &) G [0, 1]^. Then 

Cl''^'{u, v) < C{u, v) < C^/'\u, v), {u, v) G [0, 1]2, (3) 
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where 

(jOjb,8 _ niin(M, v,d + {u — a)^ + {v — b)^) 
and (j'^^'^ — niax(0, w + u — 1,0 — (a — m)^ — (6 — f )^) 

are copulas satisfying C^''''^(a, 6) = C'^^^'^ {a,b) = 9. 

Remark 1. A careful examination of the proof or Theorem 3.2.3 in [5] reveals that ([3]) also holds 
if C is a quasi-copula satisfying C(a, h) = 9. 

To close this section, we recall a well-known fact on distribution functions. Given a one- 
dimensional distribution function F{x) we define its generalized inverse by 

F-^{u) ^ inf{x e R : F{x) >u}, ue (0, 1], 

with the convention inf = +oo. If the couple {X, Y) has copula C then {X, Y) has the same law 
as {F^^{U),Fy^{V)), where ([/, V) are random variables with distribution function C. 

3 Constrained Frechet bounds 

Let S* be a compact subset of [0,1]^ and Q be a quasi-copula. We denote by Cs the set of all 
copulas C" such that C'{a, b) — Q{a, b) for all (a, b) e S, and by Qs the set of all quasi-copulas Q' 
such that Q'{a, b) — Q{a, b) for all (a, b) £ S. Define 

A^''^{u,v) := mm{u,v, min {Q(a, 6) + (u - a)+ + (u - fe)+}) (4) 

(a,b)eS 

B^'^fu, v) := max(0, u + v~l, max {g(a, b) - (a - u)+ - (b - v)+}) (5) 

The following theorem establishes that A^''^ and i?'^''? are best-possible bounds of the set Qs- 
This means that they are also bounds of the set Cs, but not in general best possible. The second 
part of the theorem gives a sufficient condition under which A^''^ or B^''^ is a copula, and therefore 
a best possible bound of Cs- As a by-product of the second part, we obtain an example of copula 
which coincides with a given quasi-copula on a given increasing or decreasing set. 

Theorem 1. 

are quasi-copulas satisfying 

B'^'^iu, v) < Q'{u, v) < A'^'^^iu, v) V(w, v) e [0, 1]^. 
for every Q' G Qs and 

A^^'^ia, b) = B^'^ia, b) = Q{a, b) (6) 

for all (a, b) G S- 

ii. If the set S is increasing then B^''^ is a copula; if the set S is decreasing then A^''^ is a 
copula- 

The proof of this theorem can be found in the Appendix. 

Example 1. This example, similar to example 2.1 in [10\ shows that if S is increasing, may 
not always be a copula. Let 5 = {(5,5) i(f'|)} and Q ^ W- Then A^'^ {^,^) = 0, and 
AS.Q (l I) ^ A^ Q (i, I) = A^ Q (|, i) = i, so that the A^^Q-volume of the rectangle [i, |]' is 
equal to — i. Similarly, if 5 is decreasing, B^''^ is not always a copula. 
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Let p : Q — > M be a mapping, continuous with respect to pointwise convergence of copulas 
and nondecreasing with respect to the concordance order on Q. We are interested in computing 
pointwise best possible bounds of the sets C := {C € C : p{C) = r} and QJ" :~ {Q e Q : p{Q) — r}. 
We denote 

yl''(u,t;) max{C(u,w)|C e C''} and u) := min{C(u, w)|C G C''} 

A''{u,v):^inscK{Q{u,v)\QeQr} and &' {u,v) := inm{Q{u,v)\Q e QT} 

for {u,v) e [0,1]2. 

For (a, 6) G [0, 1]^ and 6 G lafi ■= [W{a,h),M{a,b)], we define 

p+{aAe) ■■^piC^u'^'), P-{aAO) -.^piCl''''). 

For fixed a, b, the mappings 6 p+(a, b, 9) and 6 ^ P-ia, b, 9) are nondecreasing and continuous, 
and we define the corresponding inverse mappings by 

r I— > pl^(a, &, r) := max{0 G /q,6 : P-{a, b, 9) = r} 
r p^^(a, 5,r) := min{6' G Ia,b ■ p+ia,b,0) = r}, 

for all r such that the corresponding set over which the maximum or minimum is taken is nonempty. 
Theorem 2. Let r G [p{W), p{M)]. The hounds , and B^',B^ are given by 

A''{u,v) = A''{u,v) 

B''{u,v) = B''{u,v) 
The proof of this theorem can he found in the Appendix. 

Remark 2. This result generalizes theorems 2 and 4 in [2], which treat the cases when p is the 
Kendall's r and the Spearman's p. In these two cases, A^ and B"^ are copulas. However, in general, 
this may not be the case. Let (ai,6i) G [0,1]^, (02,62) G [0,1]^, W{ai,bi) < 9i < M(ai,&i), 
^^(02,62) < ^2 < M{a2,b2) and define 

piC) = (C(ai, 61) - 9i)+ + (C(a2, 62) - ^2)+. 

An easy computation shows that 

A"{u, v) = min(M, v, 0i + (it - ai)+ + (i; - 61)+, 6*2 + (it - 02)+ + {v ~ 62) + ), 

that is, we obtain the copula A^''^ of Equation (U) with S" = {(ai,&i), (02,62)} and Q such that 
Q{ai, hi) — 9i and Q{a2, 62) = ^2- Then, example ([Ij shows that A^ is not always a copula. 

4 Copula based pricing of multi-asset options 

We consider the problem of pricing a European-style option whose pay-off depends on the values 
of two random variables X and Y. These random variables can represent the terminal values of 
two assets (in the context of equity options) or some other risk factors which influence the value 
of the option, such as the default dates of two defaultablc bonds. 

We assume that the law of X and Y under the historical probability P is unknown, or is very 
hard to estimate, so that all information comes from the prices of traded options on these assets. 



p_\u,v,r) tfre[p{W),p-{u,v,M{u,v))] 

M{u, v) otherwise 

P+^{u,v,r) ifre[p+{u,v,W(u,v)),p(M)] 

W{u, v) otherwise 
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Under the standard assumption of absence of arbitrage opportunities in the market, the option 
pricing theory impUes that there exists a risk-neutral probabiUty Q such that the option price is 
given by the discounted expectation of its pay-ofF under Q. In practice Q is not known, and only 
some incomplete information on it can be deduced from the prices of traded options on X and Y . 

We assume that these traded options include single-asset options allowing to reconstruct the 
cumulative distribution functions Fx and Fy of X and Y . For example, if X is the price of an 
asset at time T and call options on this asset with prices Px{K) '■= E''^[e~'^'^ {X — /^)^], are 
available (where r is the interest rate and K is the strike price), the distribution function can be 
reconstructed as Fx{K) = 1 — e''^ ^^ok^^ • Similarly, if X is the default date of a defaultable bond, 
the distribution function may be reconstructed from the prices of credit default swaps on this bond 
with different maturities. 

Let the discounted pay-off function of a two- asset option be denoted by f{x,y). Its price then 
becomes a function of the copula C oi X and Y: 



7riC)=E%fiX,Y)]^ / fix,y)dCiFxix),FY{y)) 
Jo Jo 

J{F^\u),Fy\v))dC{u,v). (9) 



Jo 

It is known [71 [12] that for every 2-increasing function / such that the integral in ([9]) exists, the 
mapping C ^ i^{C) is nondecreasing with respect to the concordance order of copulas. Therefore, 
if the pay-off function / is 2-increasing, and if we know that the copula C oi X and Y satisfies 
B ^ C ^ A for two copulas A and the option price satisfies 7r(i?) < 7r(C) < 7r(^). For example, 
if no additional information on the joint law of X and Y is available, the standard Frechet bounds 
lead to 

f{F^\u),Fy\v))dW{u,v)<n{C)< f f' f{F^\u),Fy\v))dM{u,v). 
Jo Jo Jo 

Since the support of dM is the diagonal v — u and that of dW is the diagonal v = 1 — u, these 
bounds are further simplified to 



[' f{Fx\l-u),Fy\u))du<TT{C)< f 
Jo Jo 



iu))du<7riC)< / f{Fx\u),Fy\u))du. 
Jo 

However, if A and B are quasi-copulas, this method no longer applies because the integral in 
([9]) may not be well defined. The following result provides an alternative representation for 7r(C) 
which can be used for quasi-copulas, and establishes other useful properties of this mapping. We 
recall [SJ Section 4.5] that for a 2-increasing function / on [0, oo)^ which is left-continuous in both 
arguments, there exists a unique positive measure /i on [0, oo)^ such that 

fi([xi,X2) X [yi,2/2)) = f{xi,yi) + f{x2,y2) - /(a^i, J/2) - f{x2,yi)- (10) 

Proposition 2. Assume that f is 2-increasing, left- continuous in each of its arguments, and let 
the marginal laws of X and Y satisfy 

E[\f{x,o)\ + \f{o,x)\ + |/(r,o)| + 1/(0, r)| + \f{x,x)\ + |/(r,r)|] < oo. 

Then, E[\f(X,Y)\] < 00 and the mapping C H- 7r(C) is well-defined for all C, continuous with 
respect to pointwise convergence of copulas and satisfies 

7t{C) = - /(O, 0) + E[f{X, 0)] + E[f{0, Y)] 

i 

fi{dx X dj/)(l - Fxix) - Friy) + C{Fx{x),Fy{Y))), (11) 
where is the positive measure on [0, cx))^ induced by f . 



oo poo 
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fi{dxxdy)C{Fx{x),FY{Y)) 



where C is the survival copula defined by 



Ciu,v) 



u + v-l + C{l 



M, 1 — w), 



and Fx and Fy are survival functions of X and Y. 

Table[T]gives several examples of 2-asset options whose pay-offs are 2-increasing (or 2-decreasing, 
meaning that — / is 2-increasing) continuous functions. These are mainly taken from [11] , For all 
these pay-ofFs, the integral with respect to fi in formula (jlip reduces to a one-dimensional integral. 
Another important example is the function /(X, Y) = XY which is also 2-increasing, which means 
that for fixed marginal distributions, the linear correlation coefficient 



is nondecreasing with respect to the concordance order of copulas. The corresponding measure fi 
is the Lebesgue measure on [0, oo)^. 

5 Application: model-free bounds on option prices 

In this section, we derive model-free bounds on the prices of two-asset options whose pay-off 
function satisfies the assumptions of Proposition [5] when extra information about the dependence 
of X and Y is given. We give four examples corresponding to different kinds of extra information 
and different option pay-offs. 

Example 2 (The case when prices of digital basket options are known). As our first example, we 
consider an application to credit risk modeling assuming that X and Y represent the times of 
default of two corporate bonds. In this context, an important problem is the pricing of the so 
called "first to default" option with pay-off at maturity T given by = ^{xay)<t or the "second 
to default" option with pay-off — l(xyY)<T- The price of each of these options is directly 
related to the value of the copula C of X and Y at the point {Fx{T), Fy{T)): 



In view of the above, we concentrate on the "second to default" options. From the prices Pk = 
E[H^J of these options with maturities Ti, . . . ,r„, one can recover the values of the copula C of 
X and Y on the increasing set {Fx{Tk), FY{Tk))k=i,...,n- Therefore, by Theorem[Tl the copula C 
of X and Y satisfies 



p{X,Y) 



E[XY] - E[X]E[Y] 
(VarXVary)5 



E'^[H^] ^ 1 - P[X >T,Y>T]^ Fx{T) + Fy{T) ~ C{Fx{T),Fy{T)) 
E^[H^] = P[X <T,Y <T]= C{Fx{T),Fy{T)). 



B{u, v) < C{u, v) < A{u, v) V(u, v) e [0, 1]^ 



with 



A{u,v) 



min(u,w, min {Pk + {u - Fx{Tk))+ + {v - FY{Tk))+}) 



B{u,v) 



niax(0, u + V ~ I, max {Pk — {Fx(Tk) — u) 



+ 



{FY{Tu)-vr)). 
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Option type 
and/(X,y) 


increasing? 


/o°° /o°° '"('^^ ^ dy)G{x, y) 


Basket option, 
{aX + I3Y - K)+ 


-h ff Q/3 > 0, - 
if < 


sgn(a/3) / GU,i^)dz. 


Call on the minimum 
{mm{X,Y) - K)+ 


+ 


G{x, x)dx 


Put on the minimum 
{K -m:in{X,Y))+ 


+ 


G{x, x)dx 


Call on the maximum 
(max(X, Y) - K)+ 




— G{x, x)dx 


Put on the maximum 
{K -mBx.{X,Y))+ 




— G{x, x)dx 


Worst-off call 

min((X-i^i)+,(F-i^2)+) 


■ 


J^G{z + Ki,z + K2)dz 


Worst-off put 

m:m{{K^-X)+,{K2-Y)+) 


+ 


j^-^-(K^'K-)G{Ki-z,K2-z)dz 


Best-off call 

max((X-i^i)+,(y-i^2)+) 




-J^G{z + Ki,z + K2)dz 


Best-off put 

max((/^i-X)+,(jr2-F)+) 




_ j.min(K„K,) ^^^^ -z,K2- Z)dz 



Table 1: Common 2-asset option pay-off functions, and the representation of integrals with respect 
to the corresponding measure fi. The plus sign indicates that the pay-off function is 2-increasing 
and the minus that it is 2-decreasing. 
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Figure 1: Improved bound on the prices of "second to default" option as function of time to 
maturity T , when the prices of options with 2 and 3 years to maturity are known and equal to the 
price in the "Gaussian copula" model with given correlation p. Left: p = 0. Right: p = —0.7. 

and the price of any 2-asset option whose pay-off function f{x,y) satisfies the assumption of 
Proposition [5] admits the bounds 

7r(S) < 7r(C) < 7r(A). 

Since, by Theorem [TJ i? is a copula, the lower bound is sharp, while the upper bound may not 
necessarily be sharp. 

As an illustration, we have computed the upper and lower improved bounds for the prices of 
"second to default" options with different times to maturity. We assume that the marginal laws of 
X and Y are exponential with parameters Ax = 0.2 and Ay = 0.3 respectively, and that the prices 
of the "second to default" options with 2 and 3 years to maturity are known. In this example, 
these two prices are computed assuming that X and Y have Gaussian copula with correlation 
p (the Gaussian copula is the industry standard). Figure [T] shows the prices of the "second to 
default" options as function of time to maturity for two different values of p, along with the price 
in the "Gaussian copula" model and the standard Frechet bounds (without any information about 
dependence). 

Example 3 (The case when prices of all options on the maximum of two assets are known) . The 
knowledge of prices of call or put options on the maximum or the minimum of X and Y, for all 
strikes, allows to recover (by differentiation) the values of the distribution function F{K, K) for 
K >Q, or, equivalently, the values of the copula C on the increasing set {{Fx{K), Fy{K)),K > 0). 
Therefore, similarly to the previous example, the copula C oi X and Y satisfies 

B{u, v) < C(m, v) < A{u, v) V(u, v) e [0, 1]^ 

where 

A{u, v) ^ mm{u, v, mm{F{K, K) + {u - Fx{K))+ + [v - Fy{K))+}) (12) 
B{u,v) = niax(0,u + w - l,-ms^{F{K,K) - {Fx{K) - u)+ - {Fy{K) - y)+}). (13) 

To illustrate this method, we have computed the improved upper and lower bounds for the 
spread option with pay-off at date T = 1 given by /{Xt, Yt) = {Xt —Yr — K)^ . To fix the marginal 
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Figure 2: Improved bound on the spread option price as function of strike when the prices of 
all options on maximum are known and equal to the price in the Black Scholes model with given 
correlation p. Left: p — Q. Right: p = —0.7. 



laws of X and Y , we assume that Xt — Xq eyiY){axW^ — ^) and Yt — Yq exp((Tj,iyf — -f-), where 
ax ~ 0.2, (jy = 0.3, ^0 = 10 = 100 and and are standard Brownian motions. We further 
assume that the prices of all options on the maximum of X and Y are equal to the corresponding 
prices in a model where W)j^ and are jointly Gaussian with correlation p. 

Figure [5] plots the improved bounds on the spread option price as function of the strike K 
for two different values of the correlation p, along with the Black-Scholes price and the standard 
Frechet bounds. For the numerical computation of the bounds, we have taken a discrete set of 
400 strikes in (fT2)l and (fT3|) and used numerical integration to evaluate (fTTjl . which reduces to a 
one-dimensional integral in this case. 

Example 4 (The case when a single option price is known) . Assume now that the extra information 
about the dependence structure of X and Y is the expectation of a function /q which satisfies 
the assumptions of Proposition [5J p{C) :— E'^[fQ{X^Y)] — r. In this case, the price of a 2- 
asset option whose pay-off f{x,y) satisfies the assumptions of Proposition [2] admits the bounds 
Tr{B'') < 7r(C) < Tr{A'') with A" and B'' given by Theorem [2 Although A'^ and B'' are best- 
possible bounds of the set of copulas satisfying p{C) = r, if they are not copulas themselves, the 
bounds on the option price may not be best possible. 

For the actual computation of and S'' we reduce the expressions for p+ (a, 6, 9) and p^ (a, 6, 9) 
to one-dimensional integrals using the results in [SJ section 3.2.3]: 

p+{a,b,9)^ fo{F^\u),F-^(u))du+ fo{F^\u), Fy\u + b - 9))du 
Jo Je 

na+b-9 pi 

+ / fo{F^Hu),FY\u + 9-a))du+ / fo{F^\u),Fy\u))du 

J a J a-\-h—0 



9 




Figure 3: Improved bounds on the price of the option on the maximum of two assets as function 
of strike K when the price of the zero-strike spread option is known and equal to the price in the 
Black-Scholes model with given correlation p. Left: p — 0. Right: p = —0.7. 

p^{a,b,e) 

na — 6 pa 

= fo{F^\u),Fy\l-u))du+ fo{F^\u),Fy\a + b-6-u))du 

Jo Ja-0 

+ MF^\u),Fy\l + e-u))du+ fo{F^\u),Fy\l-u))du 

Ja Jl-b+0 

As the first illustration of this approach, we have computed the improved bounds on the 
price of the call option on the maximum of two assets, with pay-off at date T = 1 given by 
/(Xt, Yt) = {max{XT, Yr) — K)'^, assuming that the price of the zero-strike spread option, with 
pay-off /o(Xt, Yt) ~ {Xt — 1t)^, is known (these options are indeed often quoted in the market). 
The marginal laws of X and Y are the same as in example [3J and we further assume that the price 
of the zero-strike spread option is equal to the corresponding price in a model where and 
are jointly Gaussian with correlation p. 

Figure [3] plots the improved bounds as function of the strike K for two different values of the 
correlation p, along with the Black-Scholes price and the standard Frechet bounds. Since we now 
have much less information on the dependence of Xt and Yt than in example |31 the improved 
bounds are not as narrow as in that example. Still, when the spread option price is close to one 
of its extreme values, such as, for example in the right graph of figure [31 where we have taken 
p = —0.7, the improved bounds lead to a considerable narrowing of the price interval. In the 
numerical example, p^ and p^ were evaluated by numerical integration, their inverses were then 
computed by bissection, and a further numerical integration was performed to evaluate the bounds. 



Example 5 (The case when the linear correlation of log- returns is known) . Very often, the option 
trader does not know the full two-dimensional distribution of X and Y under Q but has a strong 
view about the risk-neutral correlation of log-returns 

E [log X log Y]~E [log X] E [log Y] 

PO = i • 

(Var(logX) Var(logy))5 

In this case, one can obtain bounds on the prices of two-asset options in the same way as in 
Example m using the function fo{x,y) — log x logy, which is 2-increasing. Figure H] plots the 
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Figure 4: Improved bounds on the price of the zero-strike spread option when the correlation of 
log-returns p is known, as a function of p. For comparison, the price in the Black Scholes model 
with the given value of correlation is shown as BS Price. 

bounds on the price of a zero-strike spread option with pay-off /{Xt, Yt) = {Xt — Yt)~^ when the 
correlation of log- returns is known, for different correlation values. As we have already observed 
in Example 01 these bounds are most useful for extreme correlation scenarios, and yield little 
additional information when the correlation is close to 0. 
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Appendix: the proofs 

Proof of Theorem]^ First, observe that A can be obtained from i? by a simple transformation: 

A^-'^{u, v) = u- B^'^{u, 1 - u) = W^{u, v), 

where the bar notation was introduced in section [21 It is therefore sufficient to prove only the 
statements involving B^'^ . 

Part i. Let us first check that B^''^ is a quasi-copula. The boundary conditions follow from 
the Frechet bounds for Q. The fact that B^''^ is increasing in each argument is obvious, and 
the Lipschitz property follows because for a family of functions {fi)i£i which are Lipschitz with 
constant 1, we have 

max/((y) <\x ~ y\+ max/j(x) and max/i(a;) < |x - y| + max/i(?/), 

i i i i 

which implies that max^ fi{x) is Lipschitz with the same constant. By proposition[T]and the remark 
after it, C^''''^^"'^^ (w, w) < Q'{u,v) for aU {u,v) e [0,1]^, (a, 6) e S and Q' G Qs- Since B'^-Q is 
the upper bound of C'^^"'^^°'"''\u,v) over {a,b) e S, we have that B^'^{u,v) < Q'{u,v). 

Let us now check the property Take (a', b') £ S. From the Frechet lower bound for Q, we 
get: 

B^-^(a', b') = max {Q{a, b) - (a - a')+ -lb- b')+}. 

(a,b)eS 

For every (a, b) G S, using the Lipschitz property of Q and the fact that it is increasing in each 
argument, we get that 

Q{a, b)-(a~ a')+ - (6 - b')+ < Q{a', b') 
Therefore, the max is attained for (a, b) = (a', b'). 



12 



Part ii. Let S be an increasing set. By adding to this set the points (0,0) and (1,1), we may 
with no loss of generahty simphfy the definition of B^''^: 

B^'''^{u, v) := max {Q(a, 6) - (a - u)+ - (6 - v)+} 

(a,b)GS 

Given that B^''^ is a quasi-copula, we only need to prove property ([T]). 

Since B^''^ is Lipschitz continuous, for every £ > 0, one can find a finite increasing set such 
that sup(„ „)g[Q IB^^'^^lu, v) — B^'^{u, v)\ < e. Therefore, it is enough to prove property ([T]) for 
a set Sn — {(fli, ^j:)}"=ii where we suppose without loss of generality that a.i < ai+i and bi < 
for i = 1, . . . ,n — 1. 

The proof will be done by induction. For n — I, property ^ is straightforward. Assume 
that it holds for Sn and let an+i > a„, 6„+i > 6„ and Sn+i ■— Sn U {(a„+i, 6„+i)}. To simplify 
notation, we write Bn ■= B^"^''^ , Bn+i ■= B^^+^'^ and Qn+i ■= (5(a„+i, 6„+i). For convenience, 
we subdivide the domain [0, 1]^ onto four sets A, B, C and D as shown in Figure [S] 

To prove that Bn+i is 2-increasing, we must show that for every rectangle R C [0, 1]^, Vb„+i (R) > 
0. However, since Vb is additive over rectangles, it is sufficient to consider only the cases RCA, 
R C B, R C C and R <Z D. By construction, on A, the function Bn+i only depends on the 
coordinate w, and therefore, Vb„+i {R) = for every rectangle RQ A. Similarly, Vb„^^ [R) = for 
R C B because Bn+i is constant on B and Vb„+i(-R) = for i? C C because Bn+i only depends 
on the coordinate v on C. It remains to consider the case R C_ D. 

Let R — [ui,U2] X [wi,W2] ^ D. We must show that 



Vb„+i(-R) = max(i?„(Mi 


Qn+l 


- (fln+l 


- Ui) 


~ {K+1 


-Vl)) 


+ max(i3„(M2, 


V2), Qn+l 


- (fln+l 


- U2) 


- {bn+1 




— max(_B„(Mi, 


V2), Qn+l 


- (a«+i 


- Ui) 


- {bn+1 


-V2)) 


- max(i3„(M2, 


Vl), Qn+l 


~ (fln+l 


- U2) 


- {bn+1 


-vi))>0 



We consider separately three cases. 

- If Bn{ui,V2) > Qn+l " (On+l " Ml) - {bn+1 ^ V2) and S„(u2,Ul) > Qn+l " {dn+l " U2) - 

{bn+i~vi) thenVB„^i{R) > Bn{ui,vi)+Bn{u2,V2)-Bn{ui,V2)-Bn{u2,vi) > by induction 
hypothesis. 

- Assume _B„(wi,W2) < Qn+i — (^n+i — ui) — {bn+1 — ^2)- Then, by the Lipschitz property of 
Bn, necessarily _B„(u2, V2) < Qn+i ~ {o-n+i — ^^2) — {bn+i ^^2), and therefore, by the Lipschitz 
property of Bn+i, Vb^+^R) = U2 - + Bn{ui,vi) - B„(u2, vi) > 0. 

— The remaining case, when _B„(mi,W2) > Qn+i — {^n+i — ui) — {bn+i — V2) and Bn{u2,vi) < 
Qn+l — (flra+i — U2) — {bn+1 — Vl), IS treated similarly to the second one. 

□ 

Proof of Theorem\^ We give the proof for the bound A'^{u,v). Since the proof is only based on 
Proposition [1] which holds in the same form both for copulas and for quasi-copulas, A^ coincides 
with A^. The proof for lower bounds B^ and B^ is similar. 

Assume r G [p(M^), p_ {u, v, M{u, v))]. Then, since 9 ^ P-{u, v, 9) is increasing and continuous, 

p{Cj^ ^ ''■*) = r and therefore, A^(u,v) > pZ^{u,v,r). On the other hand, 

{p{Q)\Q{u,v) = 9}C [p_{u,v,9),p+{u,v,9)]. 

By definition of pZ^, for all 9 > pZ^{u,v,r), p^{u,v,9) > r and therefore for every Q £ Q such 
that Q{u,v) > pZ^{u,v,r), p{Q) > r. Therefore, A'^{u,v) < pZ^{u,v,r). 
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Figure 5: Illustration for the proof of Theorem [1] part ii. 

Assume now that r > p_{u,v,M{u,v)) and let C"" := (1 - i(;)C^^"'*^^"'''^ +wM{u,v). Then 
p(C") < r, p(C^) > r (by assumption of the theorem) and since p is continuous, there exists w G 
[0, 1] such that piC") = r. Since ^"{u^v) = M{u,v) for all w, this proves that A''{u,v) > M{u,v). 
On the other hand, clearly A'^{u,v) < M{u,v) (Frechet bound). □ 

Proof of Proposition\E Since / is 2-increasing, y/([0, x] x [0, y]) = /(a;, y)+/(0, 0) — /(a;, 0)-/(0, y) 
is increasing in x and y, and therefore 

\fix,y)\^\Vfi[0,x] X [0,y])-/(0,0) + /(a;,0) + /(0,y)| 

< 1/(0,0)1 + |/(a;,0)| + \f{0,y)\ + \Vf{[0,xf)\ + \Vf{[0,yf)\ 

< C{|/(0,0)| + |/(0,a;)| + |/(a;,0)| + \f{0,y)\ + l/(y,0)| + \f{x,x)\ + \f{y,y}\}, 

for some C > 0, which implies E[\f{X, Y)\] < oo. 

Let p{dx X dy) be the law of {X, Y). By Fubini's theorem and (ITU)) we then get 

tt{C) = E[f{X, Y)] = -/(0, 0) + E[f{X, 0)] + i?[/(0, Y)] 

pOO nOO pOO pOO 

+ / lx'<xiv'<yP'{dx' X dy')p{dx x dy) 

JQ JQ Jo Jo 

= -/(O, 0) + E[f{X, 0)] + E[f{0, Y)] 

nOO nOO 

+ / p{dx' X dy')P[X > x' ,Y > y'] 
Jo Jo 

= ^f{0,0) + E[f{X,0)] + E[f{0,Y)] 

poo poo 

+ / fiidxxdy){l-Fx{x)-FY{y) + C{Fx{x),FY{Y))) 
Jo Jo 

In the last integral, the integrand is positive and bounded from above by the function 1 — Fx{x) — 
Fviy) + iiiiii(-Fx(a;), -FV(2/)), which corresponds to the copula of complete dependence and is inte- 
grable by the first part of the proposition. Therefore, the dominated convergence theorem implies 
that 7r(C) is continuous with respect to pointwise convergence of copulas. □ 
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